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NOTES  ON  MURPHY’S  METHOD 


•0 n  *-hr  nquat -jor*  from  Murphy,  Ref.  It 

z  +  (H-i/^-i^)  z  -  (M+iwT)  z  =  0  (1) 


where  i  =\/l-|z|2,  and  H,  M,  and  T  are  real,  may  be  functions 

of  | z  | 2 ,  and  co  is  real,  we  wish  to  draw  rigorously  justified 
conclusions  about  the  nature  of  its  family  of  solutions. 


In  the  first  place,  the  origin  z  =  z  =  0  is  a  singular 
point,  a  solution.  The  conditions  for  stability  are  well 
known : 


H  >  0 

(co2-4M)  H2  >  w2  ( 2T-H) 2 


(2) 


with  K,  M,  T  evaluated  at  z  =  0. 

In  the  second  place,  it  turns  out  that  under  certain 
conditions  there  is  one  (or  more)  solution  of  the  form 

z  =  re  ,  with  r  and  v  constant,  a  periodic  solution.  We 
postpone  discussion  for  a  bit. 

In  the  third  place  we  are  interested  in  the  following 
situation,  if  it  should  occur: 

There  is  a  set  of  four  real,  smooth  functions  X. (0  ,0  ) 

X  1  2 

each  periodic  in  0  and  0  with  periods  l/w  and  l/«  ,  respec- 

1  «  1  * 

tively,  and  for  each  choice  of  0  ,  6  there  exist  two  real, 

differentiable  functions  0x(t)  and  02(t)  such  that  0^(0)  = 
and  with  z  =  z  +iz,  x,  =  z  .  x,  =  z , ,  x,  =  z  ,  x^  =  z„, 
x.(t)  =  X.(fl  (t),  0„(t))  is  a  solution  of  Eq.  1  and  the  only 

X  XI  2 

one  with  x.(0)  =  X.  (0.  .0,  .  If  this  should  be  true,  we 
say  that  the  vector  function  X. (0  ,0  )  define?  a  periodic 

JL  X  a 

surface  of  Eq.  1. 
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We  know  no  way  to  establish  the  existence  of  such  surfaces 
for  Eq.  1  as  it  stands.  If,  however,  we  are  willing  to  regard 
the  nonlinearities  and  dissipative  terms  as  small,  definite 
results  can  be  established.  We  therefore  work  with 


z  -  iwz  -  Mz  +  e[(H-i/i)  z  -  (m+icuT)  z]  =  0  (3) 

where  M  =  constant,,  and  H,  m,  T  depend  on  |z[2  only,  m(0)  =  0. 
We  shall  use  this  form  in  hunting  for  periodic  solutions 
(conical  yaw)  as  well  as  periodic  surfaces  (steady  mixed 

oscillations).  Let  z  =  pe1^,  then  Eq.  3  is  equivalent  to 
the  pair: 


p  -  pep 2  +  wp'cp  -  Mp  +  e  [  (H+pp/(l~p»2  ) )  p  -  mp  ]  =  0 
p'cp  +  2pcp  -  cop  +  e[Hp<p~coTp]  +  ep 2 p cp/ (1-p2  )  =  0 


(4) 


This  set  _  does  not  depend  on  cp  itself.  We  look  for  a  pair 
p  -  PQ,  <P  -  v  such  that  according  to  Eq.  4  (with  p  =  0), 

p  =  'cp  —  0.  The  first  is  satisfied  if 

v2  -  cov  +  (M+em(p2))  =  0  (5) 


or 


v 


±  . J jo2  -4  (M+em(p  2  )  ) 


(provided,  of  course,  that  this  is  real). 

The  second  is  satisfied  if,  simultaneously 


-  uT  =  0 


(6) 


To  investigate  the 
periodic  solution, 

P  =  x2,  cp  =  V  +  X3 


behavior  in  the  neighborhood  of  such  a 
we  shall,  temporarily  set  p  =  p  +  x  , 

and  expand  system  4. 


O 
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x,  =  x. 


2p  vxn  -  k>p  x 

0  3  0  3 


X, 


+  £ 


-Hx„  -  mx,  -  m’p  x 

2  1  0  1 


2P,V 


m’x. 


P0x3 


2J2  -  "I 

l- 

V  +  o(x2 ) 

f  -  2vxn  +  wx„ 


=  i  +  ef-HP0X3  “  H'P0PX1  +  WPoT,Xx3 

+  0(x2 )  +  ep2vx2/(l-P^) 


(7) 


with  2S2  =  ±  v/co2-4(M+em)  =  2v  -  w.  It  is  well  known  that  for 
a  system  such  as 


x  =  Ax  +  0(x2) 

where  x  is  a  vector  and  A  a  constant  matrix,  that  the 
behavior  of  the  solutions  near  x  =  0,  including  stability 
criteria,  is  characterized  by  the  roots  of  the  ’secular 
equation'  det  |.A-AE j  =  0,  at  least  when  det  | A |  ^  0.  For 
simplicity,  and  to  be  able  to  compare  directly  with  pages  26 
and  27  of  Ref.  1,  we  will  specialize  to  H  =  constant,  m  =  0. 

The  secular  equation  is  then  (we  also  neglect  ep  v/ (1-p  )) 


A3  +  2eHA2  +  [4S22+e2H2]  A  -  2eo*2poT’  =  0  (8) 


The  periodic  solution  is  stable  if  and  only  if  the  real  parts 
of  all  three  roots  of  Eq.  8  are  negative,  which  in  turn  is 

2 

true  if  and  only  if  (neglecting  the  e  term) , 


H  >  0  ,  w  Q,p  T*  <  0 
H  +  cop  qT  ’/4fi  >  0 


(9) 
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(of  which  the  first  happens  to  be  a  consequence  of  -the  second 
and  third).  We  can  even  find  approximate  values  for  the  roots 


Ax  =  €copQT'/2Q, 

X  a;  ±  2iS2  -  e(H+cop  Tf/4fi) 

2  ?  3  0 


Equations  65  and  66  of  Ref.  1  lead  to  only  two  roots: 

X  =  ecop  T'/2Q 

X2  =  -  e  (H+cop flT ’/■  4S2) 


Now  let  v 
suppose  £2^0. 


=  co/  2  +  £2 ,  v  =  co/  2  -  £2 ,  4£2 2 
2 

Make  the  transformation 


2 


CO 


4M,  and 


±<P1  i  <PZ 

z  =  r  e  +  r  e 
1  2 

•  .  1V1  ^  ^ 
z  =  lr .  v  e  +ir„i^e 
XI  2  2 

with  the  understanding  that  neither  ri  nor  rg  is  zero.  Just 
as  in  Ref.  1,  we  obtain,  setting 


e 


<p  +  w 
Y1  Y2 


<P,  ~  <P* 


'61  =  Vl  +  V2  +  6  01  (ri5r2’02')  +  0(6  ) 
P1  -  V2  +  e  ®2  (ri?r2,e2)  +  0 ( e2 ) 

Ri (ri  >v2 > e2)  +  0(e2) 


0.  = 


r  = 


v  -v . 


r  = 


R.  = 


-  R  (r  ,r  ,0  )  +  0(e  ) 

J7-J7  2  1  5  2  5  2 

1  2 


(-Hv  +coT)r  + 
l  l 


p  (-Hy  +coT)cos0 

a  a 

-  m  sin  6 


(  |  rir2^1-t;2)sin02(rit,l+r2^2COSe2) 


R 


Hx'2-coT)r2  + 


1-6 
(HVj 

4-  m  sin  0 


(Ev  -coT)cos0  ' 

X  2 


j  v1r2  (v  1-P2)sin02  (r2^2-h?1^1cos02) 

^  -  — 

1-0 


(10) 
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^  ^  o  ^  2  o 

with  6 “  =  r  +  r  +  2r  r  cos  8  ,  m,  H  and  T  functions  of  <5  . 

i.  A  1  «  A 

The  right-hand  sides  of  all  four  equations  are  independent 
of  8 ^  j  so  we  may  consider  the  last  three  alone.  and 

are  of  period  2ir  in  9  .  Let 

1  f  2?r' 

H  =  — -  \  H  d@ 

1  27T  JQ 

1  f  27r 

H2  =  2tt  Jo  H  cose  d0 

1  c 27r 

Ti  =if)o  Td9 
i  c 27r 

T  =  -r-  \  T  cose  d9 

2  27T  Jo 

We  find  that 


dfl 


<Vi 


•wTx) 


(H  v  -cjT  )  r 
2  2  2  2 


d»  -  r2 


+  (H.Vdff,)  *1 


Suppose  we  can  find  a  pair  r„„,  r  „  such  that  both  R  and  R 

c  10’  20  12 

are  zero.  We  can  then  apply  a  theorem  of  Poincare  (see 
theorem  5.2  of  Ref.  2)  to  the  effect  that,  (subject  to  cer¬ 
tain  differentiability  conditions)  if  also  the  determinant 
8R. 

of  the  matrix  is  not  zero,  there  is  an  e  >0  such  that 

or .  -  l 

ii  3  ' 

for  |e|  <  e  a  solution  does  exist  for  which  r  and  r  are 

X  1  A 

periodic  functions  of  9  ,  and  therefore  of  the  time,  9  in- 

2  2 

creasing  (if  v  >  v  )  strictly  monotonically . 
x  2 

The  z  ,  z  ,  z  and  z  calculated  from  this  r  (0  ),  r  (0  ) 
1  7  2  1  2  1  2  ?  2  2 

and  the  first  of  Eq.  10  define  just  such  a  periodic  surface 

as  we  were  looking  for. 
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Let  the  matrix  dR^/dv ^  =  A.  It  can  be  shown  that  the 

character  of  the  family  of  solutions  sufficiently  near  such 
a  periodic  surface  is  determined  by  the  roots  of 
det  | A-AE |  =  0  in  the  usual  way. 

Applying  these  results  to  the  case  discussed  on  pages  28 
to  30  of  Ref,  1,  one  obtains  identical  conclusions. 

We  see  that  the  methods  of  Ref.  1  are  thoroughly  justi¬ 
fied,  in  a  perturbation  theory  sense,  insofar  as  they  deal 
with  the  'singular  points’  of  the  'amplitude  plane'  and 
their  neighborhoods,  subject  to  a  minor  modification  when 
the  critical  point  is  on  an  amplitude  axis  (r  or  r  zero). 

1  A 

It  is  possible  to  extend  this  treatment  to  cover  the 
cases  of  strongly  nonlinear  moments,  but  small  yaw,  dis¬ 
cussed  in  Ref0  3,  It  is  even  possible. to  include  all  the 
kinematic  nonlinearities  (such  as  the  f/f  term)  in  the 
zeroth  approximation  as  well.  One  proceeds  as  follows: 

Write  the  equation 

z  -  [iw+i/f]  z  -  Mz  +  e  [Hz- (m+iiJwT)z  ]  =  0 
Then  set  z  =  sinp  cos cp,  z  =  cosp  sin<p,  and  finally 

1  A  >8 

(p  =  rp  cos  p.  One  obtains  a  pair  of  equations  for  p"  and  ip , 
which  are  essentially  those  of  Ref.  4.  We  can  search  for  a 
steady  state  solution  of  this  pair  just  as  before.  We  can 
even  take  e  =  1  and  get  an  answer.  This  procedure,  and  the 
calculation  of  stability  criteria  is  carried  out  in  Ref.  4. 

Now  define 


U  =  COS P 

f  p 

f  u 

v  =  -  2  \  M  sinp 

dp  =  2  \  Mdu 

Jo 

Ji 

h  =  p  +  ip  sm  p 

+  V 

q  =  ^sm  p+wu 

h  and  q  are  closely  related  to  C  and  0o ,  respectively,  of 

1 .  ? 

Ref.  3.  We  will  find  that  both  h  and  q  are  of  order  e,  and 
both  are  well  behaved  functions  of  u,  u,  h  and  q.  Further¬ 
more 
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u2  =  (1-u2)  (h-v)  -  (q-cuu)2  =  F(u,h,q) 

If  on  the  interval  [-1,1],  F(u,h,q)  has  two  s imple  roots 
separated  .by  an  open  interval  on  which  P  >  0,  then  if  h  and 
q  were  constant,  u.  would  have  a  periodic  solution  P(t,h,q) 
of  period  l/v(h,q).  Use  that  one  for  which  u(0)  =  u  =  least 

of  the  pair  of  roots. 

Let  p(0,h,q)  =  P(i>  (h ,  q)  0  ,h ,  q) .  Now  set  u  =  p(0,h,q), 

•  9 

gU  =  y(h,q)  p(0,h,q). 

h  and  q  are  now  functions  of  0,  h  and  q,  and  0  may  be 
calculated  from 


v  £>e  =  u  =  0  p0  +  h  ph  +  q  pq 


0  =  ^(h,q)  - 


hPh+qpa 


It  can  be  shown  that  this  expression  is  well  behaved  even  at 
0=0  and  l/2  where  u  =  vpn  =0.  We  have  in  fact  obtained 

V 

a  system 

0  =  y(h,q)  +  e  ©  ( 0 ,h , q) 

h  =  e  H( 0 ,h, q)  (11) 

q  =  e  Q(0,h,q) 


with  ©,  H  and  Q  periodic  (of  period  1)  in  0.  Theorem  5.2  of 
Ref.  2  may  again  be  applied  with  entirely  similar  results. 

A  detailed  discussion  of  the  more  general  version  of  the 
problem  is  given  in  Ref.  5.  Because  it  is  possible  to  re¬ 
duce  the  problem  to  third  order  (the  last  three  equations  of 
10  or  the  set  11),  the  theory  of  periodic  surfaces  has  not, 
in  fact  been  needed,  but  only  that  of  periodic  solutions. 

On  the  other  hand,  if  the  right  side  of  Eq.  11  or  3  is  a 
constant  G  /  0,  the  full  theory  of  periodic  surfaces  is 
required  to  deal  with  the  question  of  steady  mixed  modes. 

A  discussion  of  some  simple  cases  with  G  /  0  is  contained  in 
Ref.  6. 
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